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We analyze the symmetry and the nodal structure of the superconducting order parameter in a 
cubic ferromagnet, such as ZrZn2. We demonstrate how the order parameter symmetry evolves when 
the electromagnetic interaction of the conduction electrons with the internal magnetic induction and 
the spin-orbit coupling are taken into account. These interactions break the cubic symmetry and 
lift the degeneracy of the order parameter. It is shown that the order parameter which appears 
immediately below the critical temperature has two components, and its symmetry is described 
by co-representations of the magnetic point groups. This allows us to make predictions about the 
location of the gap nodes. 
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I. INTRODUCTION 

A metallic ferromagnet is characterized by the fact that 
its electronic energy bands are split by the exchange in- 
teraction between the electrons so that the spin-up bands 
have different energies from the spin-down bands. This 
has important consequences for the symmetry and the 
gap structure of possible superconducting states. In this 
article, we study the symmetry properties of the super- 
conducting order parameter in a cubic ferromagnetic su- 
perconductor, such as ZrZn2, in the limit of weak spin- 
orbit coupling, thus complementing an earlier study by 
the authors |l| carried out in the strong spin-orbit cou- 
pling limit. 

The formation of spin-singlet Cooper pairs in a ferro- 
magnet is strongly inhibited because electrons with op- 
posite momenta and spin have energies differing by the 
exchange splitting of the energy bands [^) . Therefore we 
consider here only the case of spin-triplet pairing. In 
triplet superconductivity, the order parameter has three 
components: Ajj corresponding to the pairing of elec- 
trons in the spin-up band, A^ corresponding to the pair- 
ing of electrons in the spin-down band, and A^ corre- 
sponding to the pairing of one spin-up and one spin-down 
electron. The Ajj component is expected to be very 
small for the same reasons that inhibit the possibility of 
singlet superconductivity in a ferromagnet, and thus we 
generally neglect it. In the case of zero spin-orbit cou- 
pling, there is no coupling between the three components 
of the order parameter and thus, according to the Lan- 
dau theory of second order phase transitions, only one of 
them will become nonzero immediately below the super- 
conducting transition temperature. The turning on of a 
weak spin-orbit coupling has two effects: (i) there will 
be changes to each of the three components of the order 
parameter resulting from the lowering of the symmetry 
by the presence of the spin-orbit interaction, and (ii) the 
three components of the order parameter will be mixed 
together by the presence of the spin-orbit interaction. It 
will turn out that for the ferromagnetic magnetization di- 
rected along any high-symmetry axis, and for all possible 
symmetries of the superconducting gap function, at least 



one of the dominant components Aff and A^i has either 
line nodes or point nodes in the momentum space. These 
zeros will become deep minima in the energy gap in the 
presence of the component Ajj. The bulk of this paper 
is devoted to a detailed demonstration of these results. 

Recent discoveries of coexistence of superconductivity 
with itinerant ferromagnetism in ZrZn2 [||, and UGe2 
jlj have renewed interest to the old problem of the in- 
terplay between the two phenomena. These materi- 
als exhibit a number of peculiar properties. First, in 
contrast to all previously known examples of ferromag- 
netic superconductors, such as ternary rare-earth com- 
pounds, ruthenocuprates, etc, the same band electrons 
(d-electrons in ZrZn2, or /-electrons in UGe2) are re- 
sponsible for both the superconductivity and the ferro- 
magnetism. Second, the superconductivity occurs only 
in the ferromagnetic phase. While the exchange splitting 
of the Fermi surfaces suppresses singlet Cooper pairing, 
it was shown that the exchange by spin fluctuations can 
lead to a triplet pairing both in the paramagnetic and the 
ferromagnetic phases J^] , or to the enhancement of the su- 
perconducting critical temperature T c on the ferromag- 
netic side Q . A prominent feature of the phase diagram 
of ZrZn2 is that T c grows as pressure moves away from 
the ferromagnetic quantum critical point, which can be 
explained by the exchange-type interaction of the mag- 
netic moments of the Cooper pairs with the magnetiza- 
tion density . 

Even though the microscopic mechanism of pairing is 
not completely understood, one can use symmetry anal- 
ysis to identify the possible order parameters and de- 
termine the structure of the superconducting gap. The 
symmetry group Q of the system in the normal state is 
defined as a group of transformations which leave the sys- 
tem Hamiltonian H$ invariant. If the spin-orbit coupling 
is sufficiently strong, Q contains the operations which af- 
fect both the coordinate and the spin degrees of freedom. 
In non-magnetic superconductors, time reversal symme- 
try K is not broken, and Q = S X K X where S 
is the space group of the crystal, and U(l) is the gauge 
group ||. In magnetic superconductors, time reversal 
symmetry is broken, and Q — Sm x U{1), where Sm is 
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the magnetic space group whose elements leave both the 
microscopic charge density and the magnetization den- 
sity M invariant ^j. For example, if there is a crystal 
point group operation R which transforms M to —M, 
then the combined operation KR will be an element of 
Sm, because time reversal restores the original M not 
affecting the lattice symmetry. The combined operation 
KR is anti-linear and anti-unitary, which brings about a 
number of novel features in the symmetry analysis com- 
pared to the non-magnetic case. The symmetry prop- 
erties of the superconducting state in ZrZn 2 assuming 
strong spin-orbit coupling have been studied in Ref. M 
(see also Refs. [[0], and where various aspects 
of the theory of ferromagnetic superconductors have been 
considered). However, a rather weak magnetic anisotropy 
in ZrZn 2 P] points out that the spin-orbit coupling might 
be small, which requires a modification of the analysis of 
Ref. [jjj. A peculiar feature of ferromagnetic supercon- 
ductors, which was first emphasized by Ginzburg [ fL3|| , is 
that the internal magnetic induction in the normal state 
is always nonzero. This means that the orbital motion 
of electrons and therefore the symmetry of the supercon- 
ducting order parameter will be affected by the ferro- 
magnetic magnetization even in the absence of spin-orbit 
coupling. Another consequence is that the system under- 
goes the superconducting phase transition into a mixed 
state, even in the absence of an external field. 

The paper is organized as follows. In Sec. II, the 
normal-state symmetry groups are derived assuming that 
spin-orbit coupling is neglibly small, focusing on the cu- 
bic crystal symmetry relevant for ZrZii2. In Sec. HI, 



particle Hamiltonian 



the effect of the electromagnetic interaction on the sym- 
metry of the spin-triplet order parameters is analyzed, 
and the predictions are made about the location of gap 
zeros. The lattice periodicity is taken into account prop- 
erly, which allows us to list all possible gap nodes, in- 
cluding those at the surface of the first Brillouin zone. In 
Sec. [V, the evolution of the order parameter symmetry 



in the presence of spin-orbit coupling is studied, and it is 
shown how the order parameter is induced on both sheets 
of the Fermi surface. In Sec. [v| the Ginzburg-Landau 
free energy functionals are derived for different magnetic 
symmetries. Sec. VI concludes with a discussion of our 



results and their implications for the experiment. 



II. DERIVATION OF THE SYMMETRY GROUP 
AT ZERO SPIN-ORBIT COUPLING 



We consider the case of cubic symmetry appropriate for 
ZrZn 2 , which has the cubic laves phase structure. Also, 
we consider a single spin-degenerate electron band which 
is split by an exchange field in the ferromagnetic state. 
The symmetry of the normal (i.e. non-superconducting) 
state will be analyzed in terms of the effective single- 
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Here e is the absolute value of the electron charge, U(r) is 
the periodic crystal lattice potential, er = (<Ji, 02, 03) are 
Pauli matrices, and h ex (r) is the exchange field. The 
magnetic induction inside the ferromagnet (in the as- 
sumed long cylinder geometry) is B = curl A(r) = AttM, 
and g is the Lande g-factor for electrons, which deter- 
mines the Zeeman splitting. In the case of a collinear 
ferromagnet, which is assumed here, h ex (r) = h n f(r), 
where f(r) has the same periodicity as U(r), and Hq is 
the exchange field direction, which does not vary in the 
crystallographic unit cell. We also assume that B is uni- 
form and there is no external magnetic field (otherwise 
B = H ext + 4ttM), so that the vector potential can, 
for example, be written as A(r) = [B X r]/2. In prin- 
ciple, the magnetic induction varies both in magnitude 
and direction in the crystallographic unit cell, and our 
B is the unit cell average of the microscopic magnetic 
induction. If the variation of the magnetic induction in 
the unit cell were taken into account it would change the 
symmetry analysis given below. However, since B (ap- 
proximately 400 Gauss at zero pressure) is much smaller 
than the exchange field in ZrZii2, and since a spatial av- 
erage of magnetic induction over the unit cell is usually 
assumed to be appropriate in the calculation of the ef- 
fects of the magnetic induction on the orbital motion of 
the electrons, the approximation of a uniform B is suffi- 
cient. The exchange field h ex , the magnetization density 
M, and the magnetic induction B all have a common 
direction. 

The spin-orbit coupling in not included in Eq. (Q). 
It should be noted that when we refer to spin-orbit cou- 
pling in this article, we mean the single-particle spin-orbit 
coupling which is shown explicitly in Eq. (p5|) below. In 
principle, the microscopic magnetic dipole-dipole inter- 
action that gives rise to the internal magnetic induction 
B = AnM couples the spin and orbital motions, but be- 
cause we assume a uniform B this does not affect our 
symmetry analysis. Even in the absence of the spin-orbit 
coupling ( p5| ) , there is an effect of the ferromagnetic mag- 
netization density on the orbital motion of the electrons, 
which we refer to as the electromagnetic interaction. This 
means that the symmetry and the free energy of the su- 
perconducting state will depend on the direction of M . 

At zero spin-orbit coupling, the symmetry operations 
act independently in the real (orbital) space and the spin 
space, so that the full symmetry group of Hq is a direct 
product 



Q — G or b X G S pi n X U(l), 



(2) 



where U{\) is the gauge group composed of phase ro- 
tations <j>-0 Q ,(r)<I >_1 = e l ^tfj a (r). For our purposes, the 
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translations are not important, so that G or b contains or- 
bital rotations R or b'- 



Rorbi> a {r)Rj b = ip a (R or \r), 



and inversion: 



iMr)!- 1 =i> a (-r). 



(3) 



(4) 



Also, the effects of time-reversal symmetry are included 
in G or b- Below we shall use the notation Ckn for the 
rotations by an angle 27r/fc about the axis h in orbital 
space. The group G sp i n contains spin rotations R sp in- 

R spin Mr)R7X = [D {im {Rspin)] a ^0{r), (5) 

where DW>(R) is the spinor (j — 1/2) representation 
of rotations: for a rotation R by an angle 9 around n, 
£>(i/2)(_r) = exp[-i(9/2)(a ■ n)]. It is convenient to in- 
troduce an orthogonal basis of unit vectors ei, 6=2,(23 in 
the spin space, such that £3 \\ B. We shall use the no- 
tation Cj! n for the rotations by an angle 2-7r/fc about the 
axis n in spin space. 

A standard representation for the time reversal opera- 
tor K is K — C| e2 ifo) where Kq is the complex conjuga- 
tion operator associated with the representation {r, s z } 
[|14]| , The anti- unitary operator Kq is defined more ex- 
plicitly by the equation 



K [cip a (r)]K = c*ip a (r), 



(6) 



where c is an arbitrary c- number. In the momentum rep- 
resentation, Kq also reverses the sign of k. It should 
be noted that, in the decomposition K = C% e2 K , while 
C| e2 is an operator in spin space only, Kq operates in 
both spin and orbital space [as indicated, for example, 
by the result KqS v Kq — —s y , where s y = {h/2)a 2 is the 
y-component of the electron spin operator] . Nevertheless, 
in discussing the symmetry properties of the Hamiltonian 
Hq given by Eq. ([!]) for the case where the common di- 
rection of h ex and B is along e.3 (so that the Hamiltonian 
does not contain 02), it is useful to consider Kq together 
with the symmetry operations in orbital space. 

If R orb leaves the periodic potential U(r) and the ex- 
change field h ex (r) invariant, then the transform of the 
Hamiltonian Hq, namely RorbHoR^, is the same as 
Hq except that the vector potential A(r) is replaced by 
A'{r) = R~S b A(R orb r). This means that the transforma- 
tion rule for the magnetic induction under R orb is B(r) = 
curlA(r) -» B'(r) = curlA'(r) = i? or 1 b B(i? orf r) = 
R~ rb B(r). Also, KqHqKq is the same as H except 
that A(r) is replaced by —A(r), so that the transforma- 
tion rule for the magnetic induction under Kq is simply 
B — > — B. Thus, if R or b leaves U(r) and h ex (r) invari- 
ant, and R orb B = —B, then KQR or b is a member of the 
symmetry group of Hq . For convenience, such combined 
symmetry elements will be included with purely orbital 
elements in the definitions of the various orbital symme- 
try groups below. 



In the non-magnetic case, i.e. at h ex — B = M = 0, 
the orbital symmetry of Hq is determined by the symme- 
try of the lattice potential U(r). Since ZrZii2 has a cubic 
Laves phase structure, G or b = Oh x Ko = O x I x Kq, 
where I = {E,I} and Ko = {E, Kq}. In addition, Hq is 
invariant under arbitrary rotations in the spin space, so 
that G spm = SU(2). 

In the ferromagnetic case, where M, h ex and B are 
all nonzero, time reversal symmetry is broken, and, as 
noted above, the symmetry group of Hq contains ele- 
ments of the form KQR or b as well as purely orbital trans- 
formations. In addition, the symmetry group of Hq will 
contain operations that are purely spin-space rotations. 
More precisely, it is evident from Eqs. (|l|) and (j|) that 
Hq is invariant under the operators of the group C^ oe3 
describing the set of all spin rotations about the axis (23, 
which, as always, is taken to lie along the common direc- 
tion of M, h ex and B. Therefore, 



(7) 



This spin-space symmetry group will be combined with a 
number of orbital symmetry groups to describe a number 
of different cases corresponding to different orientations 
for the ferromagnetic magnetization density. The differ- 
ent cases will be called Case A, Case B, ... Case E. The 
appropriate symmetry groups will be described immedi- 
ately, and the order-parameter symmetries for each of the 



cases will be described later in section [II. 

Case A. The orbital symmetry of the system is de- 
termined by the electromagnetic interaction of the con- 
duction electrons with the induction B via the vector 
potential A. If this interaction can be neglected, which 
amounts to setting e — > in Eq. (|Tj), then the Hamilto- 
nian is real, so that the orbital symmetry is independent 
of M and is described by the cubic group O^, i.e. 



G orb = O/j X K = O X I X K 



(8) 



In this case, which might be appropriate for a neutral 
Fermi system, such as the liquid 3 He in magnetic field, 
or the "cold" atomic gases, only the spin symmetry is in- 
fluenced by the presence of ferromagnetic magnetization. 

In ferromagnet metals, the electromagnetic interaction 
is always present, and the presence of the magnetization 
affects the orbital symmetry even in the absence of spin- 
orbit coupling. The structure of the orbital group de- 
pends on the direction of magnetization density. In ZrZn2 
the magnetic anisotropy is sufficient weak that it should 
be possible to align the magnetization density along an 
arbitrary direction in the crystal by applying an external 
magnetic field along that direction. We now consider a 
number of possible orientations. 

Case B. If the magnetization density lies along the 
[001] direction, the orbital symmetry group is 

G or b — D 4 (C4) X I = {E, C± z , C22, C iz , 

KqC2 Xi KqC 2v , K C 2a , K Q C 2b } x I, (9) 
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where a — (x + y)/y/2, and b = (x — y)/y/2. Here we 
use a standard notation for the magnetic group G(H) 
p5[ , where the subgroup H in parentheses (the unitary 
subgroup) includes all elements of G which are not multi- 
plied by the anti-unitary operation Kq. A useful observa- 
tion is that any magnetic group G(H) can be expressed 
in terms of left cosets with respect to the unitary sub- 
group H: G(H) = H + AH, where all elements of the 
coset AH are anti-unitary. The choice of the anti-unitary 
group element A is arbitrary and does not affect the final 
results, but once chosen it remains fixed. For the group 
D^Ci), we choose A — KqC 2x . 

Case C. When the magnetization density lies along 
the [111] direction, the orbital symmetry group is 



G orb = T) 3 (C 3 )xI = {E,C 3e ,C. 

K C 2b ,K C 2 v,K Q Cw}xl, (10) 



3e J 



where e = (x + y + £)/\/3, b = C 3t b = (y — z)j y/2, and 
t>" = C 3e b = (z — x)/i/2. For this magnetic group, we 
choose A = KCih- 

Case D. When the magnetization density lies along 
the [110] direction, the orbital symmetry group is 



D 2 (C 2 ) x I = {E, C 2a , K Q C 2b , K C 2z } x I. 



(11) 



In this case, we also choose A = KC 2 b- 

Case E. For the magnetization along a general direc- 
tion, the orbital symmetry group is 



G or b — Cj — Ci x I, 



(12) 



where Ci consists of the unity operation E. This group 
is trivial and does not contain anti-unitary elements. 

In the next section, we study the symmetry properties 
of the superconducting order parameter at M using 
G spm from Eq. (jjj), and G orb from Eqs. (gggg^Q. 
The microscopic origins of the ferromagnctism and the 
superconductivity are not important for the symmetry 
analysis. 



III. SUPERCONDUCTING ORDER 
PARAMETER AT ZERO SPIN-ORBIT 
COUPLING 

In ZrZn2, the exchange band splitting is E ex ~ 
5mRy ~ 800K Ipjf , which greatly exceeds the supercon- 
ducting critical temperature T c ~ 0.2K. In this condi- 
tions, the usual Chandrasekhar-Clogston arguments ||] 
make any pairing of electrons with opposite spins, in 
particular in the singlet channel, strongly suppressed. 
The general form of a spin-triplet superconducting or- 
der parameter is A a p(k,r) — (i(ra 2 ) a pd(k,r) S. It 
is convenient to use the following representation: d(k) = 
e + d- (k) + e_d+(k) + e 3 d 3 (k) , where e± = (ei±ie 2 ) / ^/~2 
and d± — (d\ ± id 2 )j\p2. 



According to the Landau theory of phase transitions, 
the spin vector d, which appears at the critical tempera- 
ture T c , should correspond to an irreducible representa- 
tion of the normal state symmetry group Q . The easiest 
way to obtain the transformation properties of the order 
parameter under the operations from Q , i.e. the orbital 
and the spin rotations, and also the operation Kq, is to 
use the mean-field expression for the pairing Hamilto- 
nian: 



H, 



fc a,/3=l,l' 



[A Q/3 (fc)4 c 



-k,0 



H.c. 



(13) 



Here Aft = — \j2d-, which corresponds to a gap on the 
spin-up Fermi surface; An = \/2d+, which corresponds 
to a gap on the spin-down Fermi surface; and A-j-j = 
Aj.| = c?3, which corresponds to a pairing of a spin-up 
electron with a spin-down electron. Because of the Pauli 
principle, d(-fc) = -d(fe). From Eqs. |[|), (|), (|), and 
(pf), we obtain 



R orb d a (k) = d a (Rj b k) 
R spi „d a {k) = [i) (1) (]? S p m )] a( 9^(fc) 
K d a (k) =d* a (-k), 



(14) 



where a = ±, 3, and D^(R) is the vector (j = 1) repre- 
sentation of rotations. 

Since Q is a direct product of the independent orbital 
and spin symmetry groups (^), the basis functions of 
the irreducible representations of Q are given by prod- 
ucts of the basis functions of G or f, and G sp i n . An im- 
portant point here is that, because of the presence of 
the anti-unitary operations KqR 

orb in Gorb, the symme- 
try analysis should be modified. The order parameter 
should transform according to one of the irreducible co- 
representations of G or b = G(H), which can be derived 
from the irreducible representations of the unitary sub- 
group H|. 

At M = 0, G sp in — SU(2), and d transforms ac- 
cording to the three-dimensional vector representation of 
SU(2), whose basis functions are e± and e 3 . All three 
spin components d± and d 3 have the same critical tem- 
perature. At M 7^ 0, the spin symmetry is reduced to 
Gspin — Cooe 3 [see Eq. (g)], and the vector representa- 
tion is split into three one-dimensional representations of 
the group C^ e3 • The spin components d± and d 3 have 
different critical temperatures, and we assume that the 
maximum T c is achieved for d_. Thus, the order param- 
eter can be represented as an expansion 



"r 



d r (fe, r) = ie + Vi{r)fr,i(k). 



(15) 



Here fv,i{k) are the odd basis functions of a nr- 
dimensional irreducible co-representation T of G or b (the 
parity of the spin-triplet order parameter is fixed, and 
the inversion operation can be omitted from G or b)- The 
action of the orbital symmetry elements on the functions 



5 



fr,i{k) in the momentum space is denned as follows: un- 
der the crystal rotations, R or \ } f{k) — f(R„ rb k), under 
the combined operations, K R or bf(k) = f*(—R~ rb k). 
The expansion coefficients r]i(r) play the role of the or- 
der parameter components, which enter the Ginzburg- 
Landau free energy functional. The factor i on the right- 
hand side of Eq. ([lB]) is introduced so that, as we shall 
see in Sec. IV, the anti- unitary combined operations KR 



are equivalent to complex conjugation when acting on rji. 

The physical meaning of Eq. (|l^) is that the order pa- 
rameter appears only on the spin-up sheet of the Fermi 
surface, while the spin-down sheet remains normal (for 
the order parameter on the spin-down sheet, one would 
have d+ ^ 0, i.e. dae_). It should be mentioned here 
that the band structure of ZrZn2 is quite complex [ 1G 17[ . 
but we neglect such complication here and assume that 
there are only two exchange-split bands. This assump- 
tion should not affect the essence of our results. In con- 
trast to the strong spin-orbit coupling case considered in 
Ref. jij, the interband interactions c^,|C^ fc |Cfc<xc_fc/ ^, 
which could induce the order parameters of the same 
symmetry on both sheets of the Fermi surface, are absent 
due to the spin conservation. The critical temperature 
for the order parameter d^, which describes the Cooper 
pairing of electrons with opposite spins, is expected to 
be much smaller than those for d± , because of the large 
value of the exchange splitting in ZrZn2, mentioned in 
the beginning of this Section. For the same reason, we 
also neglect the possibility of a superconducting state 



with a non-zero momentum, i.e. with (c 



( L arkin- O vchinnikov- Fulde-Ferrellstate) 

When the symmetry is described by one of the mag- 
netic point groups (||), (|io|), (|TT|), or ([l2|), G or b has only 
one-dimensional co-representations (see below), therefore 
Eq. ( |l5| ) reduces to the form 



d r (k,r) = ie+r)(r)f r (k). 



(16) 



Thus, the order parameter has one component, and the 
Ginzburg-Landau functional has the same form as for 
the conventional s-wave pairing. This means that the 
phase transition from the normal ferromagnetic state to 
the superconducting state occurs into the usual mixed 
state with a lattice of the Abrikosov vortices. However, 
in contrast to the s-wave case, the orbital symmetry is 
non-trivial, in particular, there are zeros in the spectrum 
of elementary excitations where fr(k) = 0. Below we 
examine the order parameter symmetry for different cases 
and determine the positions of the gap zeros dictated by 
the magnetic symmetry. 



TABLE I: The examples of the odd basis functions for the 
irreducible representations of the point group O from Ref. 
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r 


fr(k) 


Ai 


kxkyk z (k x k-y)(M"y ^z)i}^z ' &x ) 


A 2 


kx kykz 


E 


kx ky kz ( k x ~\~ uj ky ~\~ to k z ) ? kx ky kz [k x ~\~ uj ky ~\~ ojk z ^ 


Fi 


kx 1 ky , kz 


F 2 


kxiky — k z ), k y (k z — k x ), k z (k x — ky) 



has 2 one-dimensional (A\ and A2), 1 two- dimensional 
(E), and 2 three-dimensional (F\ and F2) representa- 
tions. The examples of the basis functions are given in 
Table [j]. One-component order parameters cLa 1 (k) and 
d,A 2 (k) have line zeros at the Fermi surface, which do 
not depend on the choice of the basis functions. For 
the higher-dimensional representations, the form of the 
order parameter and its gap structure are obtained by 
minimizing the free energy in the superconducting state. 
The explicit expressions for the Ginzburg-Landau func- 
tionals and the phase diagrams for the multi-component 
order parameters can be found, e.g. in Ref. ||. 

In a charged Fermi system, where the vector potential 
created by the internal magnetization affects the single- 
electron wave functions, the cubic symmetry is reduced 
to one of the magnetic groups (|9|), ([l0|), or (p"l|), and 
the degeneracy of the two- and three-dimensional order 
parameters is lifted. Mathematically, this corresponds to 
the splitting of higher-dimensional representations of O 
into several one-dimensional co-representations. If M | 
[001] and O — > D^C/t), then it is easy to check, using 
Table H that 



Ai -» A 
A 2 ^B 

A + B 
A+ 1 E 
B 



(17) 



Fi 

F 2 



1 E + 2 E. 



We also gave here the correspondence between the 
one-dimensional representations of O and the co- 
representations of D4(C4). If M || [111] and O — > 
D 3 (C 3 ), then, using Table H, 



Case A: G or 



O x I x K 







In this case, which is relevant for the superconductivity 
in a neutral ferromagnetic Fermi system, the orbital sym- 
metry is not affected by the presence of a non-zero M . 
The order parameter is given by Eq. (15). The group O 



A x 
A 2 
E - 
Fi - 
F 2 



* A 
-> A 

1 E + 2 E 
>A + 1 E- 
► A + 1 E- 



(18) 



l E 
l E. 
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If M || [110] and O -> D 2 (C 2 ), then, using Table IFv 



A 2 

£ - 

TV 



Al 

> A - 

> A - 



B 

- B 

-A 



(19) 



B 
B. 



The physical origin of the order parameter splitting can 
be easily traced using the phcnomenological Ginzburg- 
Landau theory. For example, consider an uncharged 
Fermi liquid as above and let rj — (r/ x , r) y ,r) z ) be a three- 
component order parameter corresponding to the vector 
representation F% of the orbital group G or -b = O x I x Kq 
and corresponding to a gap function on the spin-up Fermi 
surface. Then the Ginzburg-Landau free energy describ- 
ing a homogeous phase is 

F = arj* -T7 + /3i(>7* • r]) 2 

+P2\ri-v\ 2 + M\rix\ 4 + \Vy\ i + \Vz\ i ), (20) 

where a = a(T— T Cj q), and T Cj o is the critical temperature 
at e = 0. There are a number of physically different 
states that minimize this free energy, depending on the 
values of the parameters of the fourth-order terms [ jl9| ; 
for example, one of these solutions has the form rj — 

Now, for a charged (metallic) ferromagnet, it is impor- 
tant to include the gradient terms in the free energy, so 
that the terms in the free energy of second order in the 
order parameter becomes EQ] 



F = 



a(T - T Cj0 )|i7| 2 + Ki(D i ri J )*(Di'qj) 
+K 2 (D i r h )*(D j r )j ) + # 3 (Ar?i)*(£>^i) 
+K 4 (D i r h )*(D i r H ) 

a(T - T Ct0 )\v\ 2 + ijfo* X V\B + K^D^)* {D %m ) 



-K 23 [(D lVi y(D jVj ) 
-K 4 {D m y{D m ). 



(D lV3 y(D jr]l )} 



(21) 



Here D = V + z(27r/<J> )A, <j> = nhc/e is the flux quan- 
tum, K 23 = (K 2 + K 3 )/2, and 7 = ir(K 3 - K 2 )/$o- In 
the second part of Eq. (pll), we regrouped the gradient 
terms using the identity ^Dj, Dj] = — (27ri/$o)eyfc-Bfc. 
The quantity i[rj* x 77] can be interpreted, up to a fac- 
tor, as the density of the orbital magnetization of Cooper 
pairs The second order terms given by Eq. (|2l]) are 
sufficient to calculate the critical temperature describing 
the phase transition from the normal state to the super- 
conducting mixed state. The free energy of the super- 
conducting state will depend on the direction of the flux 
lines (determined by the direction of M relative to the 
underlying crystal lattice). 

Here we consider only the case M || [001], so that 
G or b = D4(C 4 ). The critical temperature for the order 
parameter component rj z can be calculated exactly, while 
the critical temperatures for rj± = rj x ± ir\ y can be found 



TABLE II: The character table and the examples of the odd 
basis functions for the irreducible co-representations of the 
magnetic point group D^Czi). The overall phases of the basis 
functions are chosen so that KoC2xfr(k) = fr(k). Ai,2 are 
arbitrary real constants. 



r 


E 


Caz 


fr(k) 


A 


1 


1 


k z 


B 


1 


-1 


fc z [Ai(fc !/ + ik x ) 2 + A 2 (/c !/ — ik x ) 2 } 


l E 


1 


i 


ky ~\~ ikx 


2 E 


1 


—i 


ky ikx 



using the variational approach similar to that of Ref . [£lj , 
with the result: 



Tc,+ 


— T C fi 


-t— I K\ -\ 


-K 3 + 


Ki 
2 


Tc,- 


= 7c, 




-K 2 + 




T c , z 


= 7c, " 


a$ 







M 



Barring accidental degeneracies, these critical temper- 
atures are all different, so that, at e ^ 0, the three- 
component order parameter is split. The difference be- 
tween the critical temperatures T Ci+ and T c ._ is propor- 
tional to 7, and is entirely due to the interaction of the 
orbital pair magnetization with B. It is easy to see, using 
Table O, that the order parameter components 77+ , ?/_ , 
and rj z correspond to the following one- dimensional co- 
representations of D^Cj): rj z ~ A, rj + ~ 1 E : rj- ~ 2 E. 
It may be that as the temperature is lowered below this 
critical temperature into the superconducting state and 
the fourth-order terms in the free energy become more 
important, there will be a second phase transition that 
into a state that does a better job of minimizing the 
fourth-order contributions to the free energy. 



Case B: G orb = D 4 (C 4 ) x I (M \\ [001]) 

The order parameter is given by Eq. (|l6|), and the 
irreducible co-representations are listed in Table ||. We 
see that the order parameters d,A(k) and de(fc) vanish 
on the line k z = at the Fermi surface, while di E (k) 
and d,2 E (k) vanish at the points k z — k y — [note that 
here the label T refers to the orbital symmetry, whereas 
111 Ref. § we labelled the co-representations by their 
total (orbital plus spin) symmetry]. These zeros are not 
accidental in the sense that they are independent of the 
choice of the basis functions. Indeed, one of the elements 
of the magnetic point group D 4 (C 4 ) is the two- fold ro- 
tation C 2z . Therefore, 



C 2z fA,B{k) = fA,B(-k x 

= ~ fA,B{k X T ky, —k z ) = 



ky, kz) 

fA,B[k), 



(23) 



7 



so that fA,B(k x ,k y ,0) = 0. Similarly, under a four- fold 
rotation around the z axis: 

C± z f B {k) = fB{k y ,-k x ,k z ) = -f B {k), 

therefore /b(0, 0, fe z ) = 0. Also, 

Ci z f\E?E{k) = fiE, 2 E(ky- ~k x ,k z ) = ±ifi E 2 E (k), 

hence fi E ^E{0, 0, A;*) = 0. 

It also follows from Eq. ( p3| ) that /^(fe) and /b(&) go 
to zero at the surface of the Brillouin zone, i.e. at k z — 
±ir/a (a is the lattice constant), because (k x ,k yi Tr/a) 
and (k x , k yi — 7r/a) are equivalent points. In order to 
take into account the crystal periodicity leading to the 
presence of these additional gap zeros, one has to rep- 
resent the basis functions as the lattice Fourier series 
f(k) = J2 n fn e%k R " ) where summation goes over the 
sites Rn of the Bravais lattice of the crystal. The ex- 
pansion appropriate for an odd order parameter has the 
form 



f(k) = ^c„sinfc • Rn, 



(24) 



where R n are the sites of a fee cubic lattice, 
which cannot be transformed one into another by 
inversion. In the nearest-neighbor approximation, 
we choose the following set of JR„'s: {R n } — 
(a/2H(101), (101), (011), (011), (110), (110)}. Using Ta- 
ble |l|, we obtain the basis functions which have 
symmetry-imposed zeros at the surface of the Brillouin 
zone: 



TABLE III: The character table and the examples of the 
odd basis functions for the irreducible co-representations of 
the magnetic point group D3(Cs). The overall phases of the 
basis functions are chosen so that ifoC?2ii/r(fe) = fr(k). 



r 


E 






Hk) 


A 


1 


1 


1 


kx ~\~ ky -\- k z 


l E 


1 


Id 


LJ* 


O rv> x y 2 


2 E 


1 






C t\ix Ivy \^ C ISjz 



TABLE IV: The character table and the examples of the 
odd basis functions for the irreducible co-representations of 
the magnetic point group D2(C2). The overall phases of the 
basis functions are chosen so that KoC2b.fr(k) — fr(k). A is 
an arbitrary real constant. 



r 


E 




Hk) 


A 


1 


1 


kx ~h ky 


B 


1 


-1 


k z + i\(k x — k y ) 



Case C: G orb = D 3 (C 3 ) X I (M || [111]) 

The order parameter is given by Eq. (|l6|), and the 
irreducible co-representations are listed in Table 



HI. The 



order parameters di E (k) and d2 E {k) vanish at the points 
where the line k x = k y = k z crosses the Fermi surface, 
but d^(fc) does not have zeros. The zeros of diE 2 E(k) 
are imposed by symmetry, because under a three-fold 
rotation about the axis e, 



f A (k) = sin- 



k,a 



cos ■ 



cos 



k z a 



kr 



k t . 



fs(fe) = sin cos cos — — 

J v ' 2 V 2 2 

p \ k z a f . kyd . , k x a 

fi E {k) = cos — I sin ~y + 1 sin — 



+Ai 



i tv . / k x a 


fc,,et\ 


e 4 sin 1 — - — 







— e 4 sin 



, , k z ci ( , kyd , . k x ci 
PeW = cos — I sin -| i sin — 



+A 2 



4 sin 



2 



e 4 sin 



2 



Here Ai^ are arbitrary real constants. The polynomial 
expressions for the basis functions from Table [f| are re- 
covered in the limit of a "small" Fermi surface k — > 
[note that /s(fc) from Table [n] can be obtained by includ- 
ing the next-nearest-neighbors in the expansion (p4|)]. It 
should be noted that these nearest-neighbor results give 
also gap zeros not required by symmetry, e.g. /s(fc) = 
on the plane k x = k y . These "accidental" zeros will be 
removed if higher-neighbor terms are included, but if the 
nearest-neighbor terms turn out to be dominant, experi- 
ment could find indications of these accidental zeros. 



C3 e /i Ej 2 E (fc) = /i Bi 2 E (fc z , k x , k y ) = e 



_ ±27ri/3 



fiE?E(k), 



so that fi E , 2 E(k x — k y = k z ) = 0. 

We also give the expressions for the basis functions of 
the magnetic point group Da(C3) in terms of the lattice 
Fourier series in the nearest-neighbor approximation: 



f A {k) = 
h E {k) = 

f' E (k) = 
where Sf 



S+ + S+ + S+ + iAi (5f + Sz + S 3 -) 
uj*S+ +ujS+ + S+ 

ujS+ +uj*S+ + S+ 

+1X3 (UlS-~ + U!*S2 + S3) , 

sin(fc x a/2 ± k y a/2), S^ 1 — sm(k y a/2 ± 



k z a/2), S 3 = sin(fc z a/2 ± k x a/2), and are arbi- 

trary real constants. 



CaseD: G or6 = D 2 (C 2 ) x I (M || [110]) 
The order parameter is given by Eq. (|l6|), and the 



irreducible co- representations are listed in Table IV. The 
order parameter ds{k) does not have zeros, but ^^(fc) 
has the symmetry-imposed lines of zeros where the plane 
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k x = —k y crosses the Fermi surface, because under a 
two-fold rotation about the axis a, 

C^2o/A(fe) = /a(^j/i ^i! ~kz) = — Ja{~ ky, —k x: k z ) 
= fA(k), 

so that fA{k x = — fcy) = 0. 

The basis functions of the magnetic point group 
D3(C 3 ) in terms of the lattice Fourier series in the 
nearest-neighbor approximation: 



f A (k) = cos — I sin — + sin -|- 



. . 1 k X 0- kyCl 

-Ai sin h - JL - 

1 2 2 



/ B (fe) = sin — I cos — + cos , ( 
+iA 2 ainf — 



where A1.2 are arbitrary real constants. 



Case E: G or b = Ci x I 

The group Ci has single one-dimensional odd repre- 
sentation, which is realized by any odd function of k. 
Therefore, there are no symmetry-imposed gap nodes in 
this case. 



IV. SUPERCONDUCTING ORDER 
PARAMETER AT WEAK SPIN-ORBIT 
COUPLING 

Now let us turn on a weak spin-orbit coupling neglected 
in the previous discussion. We shall see that the effect of 
spin-orbit coupling is two-fold. First, it mixes together 
the order parameters on different sheets. Second, similar 
to the electromagnetic interaction studied in the previous 
sections, it reduces the symmetry of the order parameter 
and changes the gap structure on each sheet of the Fermi 
surface. 

In the presence of spin-orbit coupling, the normal state 
Hamiltonian (fil) contains an extra term: 



Hq s-o — Hq 



4m 2 c 2 



VU(r) x [p+ -A 

" c 



(25) 



Spin is no longer a good quantum number and should be 
replaced by pseudospin [p2|. In contrast to Eq. (Q), the 
symmetry group of Eq. ( p5|) cannot be represented as a 
product of independent orbital and spin groups. Instead, 
we have, neglecting the translations, 



where G S _ G consist of rotations which affect both the or- 
bital and the pseudospin degrees of freedom: 

i?^ a (r)i?- 1 = [DW i \R)\ af) il>f ) {R- 1 T\ (27) 

and also the combined operations KR, where K = 
C% e2 K , so that 



K[cip a (r)]K 1 = c*(io 2 ) a {3il>i3{r), 



(28) 



where c is an arbitrary c- number (note that K 2 = — 1). 
The transformation rules for the order parameter become 
[cf. Eq. @] 



Rd a (k) = [DV{R)] afj d {R- l k) 
Kd a {k) = -d* a {-k) = d* a {k). 



(29) 



In this shown in Rcf. jjj , the symmetry of the sys- 

tem is reduced to a magnetic point group G(H), and the 
superconducting order parameter transforms according 
to one of the one-dimensional irreducible co- representa- 
tions. Depending on the direction of the magnetization, 
G(H) = D 4 (C 4 ), D 3 (C 3 ), D 2 (C 2 ), or d (in Rcf. §, 
only the first two cases were studied). The only differ- 
ence from the previous Section is that the elements of the 
magnetic groups now act simultaneously on the orbital 
and the spin coordinates, see Eqs. ( p7| ) and (p8|), and one 
should replace KqR with KR in the definitions (|9|), jlf^), 
and (0). 

Because of the possibility of the interband pairing in- 
teractions of the form cjLc_ fe ^Ck>iC—k',l> the supercon- 
ductivity is present on both sheets of the Fermi surface. 
Instead of Eq. (|l5|), we have the following general ex- 
pression for the order parameter: 

d(k) = e+d-(k) + e_d + (fc) + e 3 d 3 (k) 

n T _ n T+ 

~ ie + ^2 V-,ifr-,i( k ) + ie~ »7+,»/r+,»( fe )-( 30 ) 



Here T_ (r+) label the irreducible co-representations of 
G(H) describing the orbital symmetry of the order pa- 
rameter at the pseudospin-up (pseudospin-down) sheets 
of the Fermi surface. The choice of these representations 
is not arbitrary because e+d- and e_d-|_ should have 
the same symmetry properties. Thus, the order param- 
eter has nr + + nr_ components: (r)+,rf-), where rj + — 
(J7+,1, ...,77+, nr+ ) and r) = (r)-,i, 7?-,„ r _ ). For the 
magnetic groups of interest to us, all co-representations 



are one-dimensional, so that nr 4 



1. As dis- 



cussed in Sec. Ill 



g = g s _ x u(i), 



(26) 



the contribution proportional to e 3 is 
small because of the large exchange band splitting, and 
is neglected in the second line of Eq. (|3(J). 

It is instructive to study the evolution of the order pa- 
rameter symmetry in the presence of spin-orbit coupling 
using the Ginzburg-Landau theory. Let us start by look- 
ing at the first of the effects mentioned in the beginning 
of this Section (i.e. the order parameter mixing), us- 
ing as an example the vector representation F\ of O and 
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assuming M || [001]. For the moment, we neglect the 
electromagnetic interaction and omit the gradient terms 
in the free energy. At zero spin-orbit coupling, the orbital 
symmetry is cubic, and T + = T_ = F±. It is convenient 
to use the following set of the basis functions of F\ : 



h{k) = 



k * +lk \ / 2 (fc) = ^-^-, /,(*) = *,, (31) 



V2 



V2 



then rj± — (r)± t i, r)±,2, ^±,3)1 and the quadratic part of 
the free energy is 



F = a+(T - T^lv+l 2 + a_(T - T c ,_)|r,_ 



(32) 



The critical temperatures T c ._ and T c + for the spin-up 
and the spin-down order parameters are different, in gen- 



> T t 



c,+) 



There are no mixed 



eral (we assume that T Ci . 
terms of the form r\\^r\~,j + c.c. in Eq. (|32|), because of 
the spin rotation symmetry U(l). Indeed, under a spin 
rotation by an angle 9 about e 3 , we have d± — > e ±l6 d±, 
which can be interpreted as an operation acting on the 
order parameter components: T]± — > e~^ % ®r\±. The mixed 
terms are not invariant under such transformations and 
therefore are not allowed. This is, of course, the same 
continuous symmetry which is responsible for the spin 
conservation. 

Now, if a weak spin-orbit coupling is turned on, we can 
treat it as a symmetry-breaking perturbation in the phe- 
nomenological Ginzburg-Landau functional. The spin ro- 
tations are no longer symmetry elements on their own, 
and, in addition to the terms on the right-hand side of 
Eq. d32|), the free energy should contain other invari- 
ants built from the components of (r; + ,?7_). The mag- 
netic group D4(C4) is generated by the rotations C^ z 
and the combined operations KC 2x - According to ( |29| ), 
C 4z d ± (k) = ±id ± (C^k), KC 2x d±(k) = d* ± (C^k). In 
terms of rj± , we have 



C 4 zV±,3 = ±i»7±,3 
KC 2x V± = V±- 



(33) 



(Note that, because of our choice of the basis functions 
and the presence of the overall factors i on the right- 
hand side of Eq. ( p9[ ) , the action of KC 2x on the order 
parameter components is equivalent to complex conjuga- 
tion.) Using Eqs. (|33|), we obtain quadratic terms which 
are invariant under all transformations from DjfCj) and 
should therefore be added to the free energy (p2|): 

3 

F s - o = F + (*+,i\v+,i\ 2 + *-,i\V-,i\ 2 ) 

i=l 



(34) 



The coefficients X±j and 7^2 are small at weak spin- or- 
bit coupling. The model of Eqs. (|32"|) and (pi|) can have a 



rich phase structure, depending on the relation between 
the "bare" critical temperatures T Cj _ and T c _ + and other 
parameters. In order to work out the whole phase di- 
agram and the structure of successive superconducting 
phases, one should include fourth-order terms in the free 
energy (^||) and (|34|), which we shall not do here. In- 
stead, we concentrate on finding the maximum critical 
temperature. 

The components (77+, 1, ^-,2), (v+,2, V+,3, and 

?7_ 3 can be considered separately. For example, the crit- 
ical temperature for (77+. 2, V-.i) is given by 




T-,i] 



4 7l 2 



(35) 



where T± i = T c ,± — A± t i/a±. Both components 77+ 2 and 
are non-zero below T c , so that superconductivity ap- 
pears simultaneously on both sheets of the Fermi surface. 
The order parameter can be obtained from Eq. (p0|): 



d(k) 



kii \ tk'x 



V2 



-rf-^i 4- ie- 



ikx 



V2 



-V+,2- 



(36) 



At weak spin-orbit coupling and T c ._ > T C) +, 77+^ is 
much smaller than 7?-,i: 77+, 2/77-, 1 oc 71. The order pa- 
rameter ([3(|) has point nodes at the poles of the Fermi 
surfaces and, according to the classification of Ref. [El, 
corresponds to the irreducible co-representation A of 
D 4 (C 4 ). Similarly, one can derive T c for the order pa- 
rameter T)—,2) and check that it corresponds to the 
co-representation B. 

The critical temperatures for r)± t 3 are T± ,3 = T c> ± — 
A±,3/a±. The corresponding order parameter d still van- 
ishes on one of the sheets of the Fermi surface, which is an 
artifact of our model, based on the representation F\ of 
O. If one includes all representations of the cubic group 
in the free energy (|3^) , then the spin-orbit coupling would 
lead to the appearance of a variety of quadratic terms 
which mix together different representations on different 
sheets, similar to Eq. (|34|). In this case, the order param- 
eter will always be present on both sheets of the Fermi 
surface, and the results of Ref. [Q will be recovered. 

Now we study how the nodal structure of the super- 
conducting order parameter on a single sheet (say, the 
pseudospin-up sheet) evolves with spin-orbit coupling. 
We consider only the case M || [001], neglect the elec- 
tromagnetic interaction, and start from the representa- 
tions Ai and F\ of the group O at zero spin-orbit cou- 
pling. The order parameters corresponding to A\ is 
d Al (k,r) = ie+i{r)f Al {k) (sec Eq. ©). The order 
parameter corresponding to F\ has the form (|l5|) with 
r = Fi, n-p = 3, and the basis functions given by Eqs. 
d3l|). The quadratic part of the Ginzburg-Landau func- 
tional is 



F = a Al (T - T Al )\H\ 2 + a Fl (T - T Fl )|r,| : 



(37) 

There are no mixed terms in Eq. (|37]) because of the 
different transformation properties of £ and r\ with re- 
spect to the elements of the cubic group. We assume 
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T Al > Tp ± , so that only £ is nonzero immediately below 
the critical temperature. From Table |[ the order pa- 
rameter d Al has six line nodes where the planes k x — 0, 
ky — - k z — - k £ — ky ; ky ~~ k z ^ tind k - — - k% cross 
the Fermi surface. However, according to Table || all 
these gap nodes, except from that on the plane k z = 0, 
are incompatible with the magnetic symmetry D^Ci). 
Let us now see how the extra nodes disappear when the 
spin-orbit coupling is taken into account. 

The spin-orbit coupling reduces the cubic symmetry 
to D^Cj), whose action on the components rj(= r}—) is 
given by Eqs. (|33|), and on £ by 



C 4z £ = -t£ 
KC 2x £, = r 



(38) 



[here we used Eq. (£9|) and the identities /a 1 (C 42 , fe) = 
f Al {k) and f* Al {~ C 2x k ) = -fAiity]. Since the compo- 
nents £ and 773 have the same transformation properties 
under all operations from D^C/i), the free energy, which 
is invariant with respect to the magnetic group, should 
contain mixed terms in addition to (BT 



Fo+liCm+r&Z), 



(39) 



where 7 is small at weak spin-orbit coupling. The critical 
temperature is changed compared to T Al : 



(40) 



and the order parameter on the pseudospin-up sheet now 
has the form 




d(k) = ie + [£f Al (k) 



■V3fF u a(k)] oc e + k z 



(41) 



This order parameter corresponds to the the co- 
representation A of 1)4(04). Thus, the only line node 
that survives the presence of the spin-orbit coupling is 
located on the plane k z = 0. However, if the spin-orbit 
coupling is weak, then the subdominant component 773 
is small, and the other five line nodes of f Al (fe) are just 
slightly filled, so that we shall have deep minima in the 
gap. At not very low temperatures, these "quasi-nodes" 
cannot be distinguished experimentally from true line 
nodes. 



V. GINZBURG-LANDAU THEORY FOR 
FERROMAGNETIC SUPERCONDUCTORS 

We have seen in the previous sections that both the 
electromagnetic interaction and the spin-orbit coupling 
break the cubic symmetry, lift the degeneracy of the order 
parameter, and change the gap structure. In addition, 
the spin-orbit coupling induces non-zero order parame- 
ters on both sheets of the Fermi surface. The symmetry 
is reduced to a magnetic group D„(C„) (n = 2,3,4), 



TABLE V: The pairs of orbital co-representations corre- 
sponding to the same symmetry of the order parameter ( fl^ ) 
on both sheets of the Fermi surface. 



G(H) 


(T+,r_) 


D 4 (C 4 ) 


(A,B), (B, A), CE, 2 E), CE^E) 


D 3 (C 3 ) 


(A, 2 E), CE,A), CE^E) 


D 2 (C 2 ) 


(A, B), (B, A) 



or Ci. All co-representations of these groups are one- 
dimensional, so that the general form of the order pa- 
rameter is given by 

d(k,r) =ie + f r _(k)j 1 -(r)+ie_f r+ (k)r 1+ (r). (42) 

The order parameter symmetry should be the same on 
both sheets, which means that (i) both components T)- 
and 77+ have the same transformation properties under 
the action of the magnetic group elements, and (ii) there 
are some restrictions as to the choice of r + and T_ , stem- 
ming from the different transformation properties of the 
spin vectors e + and e_. In Table the pairs of or- 
bital co-representations giving rise to the same symme- 
try of d are listed for all three relevant magnetic groups. 
For instance, the order parameter (^6|) corresponds to 
(r+,r_) = i^E^E). The examples of the basis func- 
tions /r ± (fe), which have only the zeros imposed by sym- 
metry, can be found in Tables |H , III, and |y|. It is easy to 
see that the order parameter always has nodes, at least 
on one of the sheets of the Fermi surface. 

The Ginzburg-Landau functional contains all possible 
uniform and gradient terms which are (i) invariant with 
respect to G(H), and (ii) gauge invariant. The uniform 
terms have the same form for all three magnetic groups: 



F, 



uniform = a+ (T - T+) \t] + | 2 + <2_ (T - T_ ) |r/_ | 2 

+l{rf+V- + rf-V+) + Fi, (43) 



where F4 is given by 

F4 = ftlTy+l 4 -)-/^! 4 

+fa\ v+ \ 2 \v4 2 + MrtW- + V-W+) (44) 

+{^\v+\ 2 + ^\v-\ 2 M+i- + i-i+)- 

The coefficients 7, f3±, (3§, /3q vanish at zero spin-orbit cou- 
pling, due to the spin rotation symmetry. 

The gradient terms are different for different magnetic 
groups. For G(H) = D 4 (C 4 ), 

F gr ad = K+lD^+l 2 + K+\D zV+ \ 2 
+K^\D^\ 2 + K^\D Z ^\ 2 
+K4(D ± r) + y(D ± r]-) + c.c.} 
+K 6 [(D z r) + )*(D z n-)+c.c.}, (45) 

where D± — (D x ,D y ). The coefficients K4 and Kq van- 
ish in the absence of spin-orbit coupling. 
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In the case of G(H) = D3(Cs), it is convenient to 
make a change of coordinates after which z is directed 
along [111]: r — > r 1 = Rr, where R is the matrix of a 
three-dimensional rotation by an angle 9 = arccos(l/V3) 
about the axis B. Omitting the primes, the gradient 
terms in the new coordinates have the same form as Eq. 

Finally, for G(H) = D2(C2), it is convenient to rotate 
the coordinates in such a way that z is directed along 
[110]: r — > r' = Rr, where R is the matrix of a three- 
dimensional rotation by an angle 8 — ir/2 about the axis 
B. In this case, the gradient terms have the following 
form: 

F g rad = K+\D xV+ \ 2 + K+\D y r,+ \ 2 + K+\D z r,+ \ 2 
+Ki\D xV ^\ 2 + K;\D yV _\ 2 + Kz\D z ^\ 2 
+ [K 4 (D xV+ )*(D x r]-) + K 5 (D vV+ y(D yV _) 
+K 6 (D z r) + )*(D zV -) + c.c.}. (46) 

Because of the choice of coordinates, M — Mz and B = 
Bz in all three cases. The coefficients K4, K5, Kq vanish 
in the absence of spin-orbit coupling. 

If the ferromagnetic magnetization is not directed 
along a high symmetry axis, then G(H) = Ci. In this 
case, the only symmetry element is the unity operation, 
and the gradient terms contain all possible real combina- 
tions of the components of D and r/± . We shall not give 
these (rather cumbersome) expressions here. 

The Ginzburg-Landau functionals listed above can be 
used for deriving the phase diagram of a cubic ferromag- 
netic superconductor, which can be quite complex. In 
particular, one cannot exclude the possibility of extra 
phase transitions in the superconducting state. 

VI. CONCLUSIONS 

We have studied the symmetry of the superconducting 
order parameter in a cubic ferromagnetic superconduc- 
tor. An experimental example is provided by ZrZn2. Be- 
cause of the anti-unitarity of time reversal, the usual sym- 
metry analysis of possible superconducting states (see 
Refs. jF9|, [ p0( | . and p|) is not applicable. In a metallic 
ferromagnet, when both the electromagnetic interaction 
and the spin-orbit coupling are present, the order param- 
eter symmetry evolves from that appropriate for the cu- 
bic group O/j to one of the magnetic point groups, which 
is studied here using the phenomcnological Ginzburg- 
Landau theory. It is shown that the order parameter 



corresponds to one of the irreducible co-representations 
of the magnetic group, and has two components, which 
describe pairing on the exchange-split sheets of the Fermi 
surface, see Eq. ([42|). It should be noted that our results 
follow from general symmetry considerations and do not 
depend on the nature of ferromagnetism in the normal 
state (itinerant vs localized moments) or the mechanism 
of superconducting pairing. 

We have determined the fc-dependence of the order 
parameter imposed by the magnetic symmetry for all 
possible directions of the ferromagnetic magnetization, 
see Table [v| The most remarkable result is that there 
should always be zeros in the energy gap, either point 
nodes or line nodes or both, at least on one of the sheets 
of the Fermi surface, when M is directed along any of 
the high symmetry axis of the cubic lattice. These nodes 
should give rise to a power-law behavior of the thermo- 
dynamic and kinetic characteristics It is expected 
that such experimental techniques as ultrasonic attenu- 
ation measurements in the superconducting state might 
be especially useful in determining the detailed structure 
of the order parameter (a discussion of this can be found 
in Ref. It should be noted that, if the electromag- 

netic and the spin-orbit interactions are weak, then the 
gap nodes appropriate for the underlying cubic symme- 
try would manifest themselves as deep minima of the gap, 
which would also have to be taken into account when an- 
alyzing the experimental data. 

The situation might be complicated by the presence 
of additional phase transitions in the superconducting 
state, which is a common feature of the systems with 
multi-component order parameters. Because of the com- 
plexity of the Ginzburg-Landau functionals derived in 
Sec. |v|, the number of possible scenarios with different 
predictions for experiment is quite large. In our view, 
it is still premature to discuss specific models, because 
of the lack of experimental data in the superconducting 
phase of ZrZn2. 
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